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Abstract 

The Breit equation with two electromagnetic form-factors is studied to obtain a potential with 
finite size corrections. This potential with proton structure effects includes apart from the stan- 
dard Coulomb term, the Darwin term, retarded potentials, spin-spin and spin-orbit interactions 
corresponding to the fine and hyperfine structures in hydrogen atom. Analytical expressions for 
the hyperfine potential with form factors and the subsequent energy levels including the proton 
structure corrections are given using the dipole form of the form factors. Numerical results are 
presented for the finite size corrections in the IS and 2S hyperfine splittings in the hydrogen atom, 
the Sternheim observable D21 and the 2S and 2P hyperfine splittings in muonic hydrogen. Finally, 
a comparison with some other existing methods in literature is presented. 

PACS numbers: 03.65.-w, 32.10.Fn, 13.40. Gp 
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I. INTRODUCTION 



The Breit equation [l, 2| is a paradigm example of how one derives coordinate potentials 
from Quantum Field Theory: an elastic scattering amplitude, expanded in powers of 
and depending on the three momentum transfer q, gets Fourier transformed into the coor- 
dinate space. The result is the potential K(r). Famous examples of this general principle 
include, among others, (i) the Casimir- Polder forces between neutral atoms as van der Waals 
forces at large distances from a two photon exchange amplitude [3j, (ii) the Uehling-Serber 
potential from one- loop vacuum polarization diagram j4|, (iii) the Feinberg-Sucher two neu- 
trino exchange force |5|, (iv) the microscopic potential of Nuclear Physics based on a, p 



and u exchanges jo] and (v) quantum corrections to the Newtonian potential |7|]. Other 
potentials are derived from pseudoscalar (Goldstone bosons) exchanges js], scalar Higgs ex- 
changes (Higgsonium) {q] and even from Finite Temperature Quantum Field Theory which 



gives temperature dependent potentials 



10| . The span of the applications of these poten- 



tials, derived via the Fourier transform of an elastic amplitude, ranges from atomic and 



molecular physics up to nuclear-particle physics and even gravitation and cosmology 

The Breit equation follows the very same principle for elastic e~fi~^, e~^e^ (positronium), 
e~p (hydrogen) and p~p (muonic hydrogen) amplitudes. The outcome is a potential which 
goes beyond the Coulomb potential and includes the correct expressions for fine structure, 
hyperfine structure and the Darwin terms. To appreciate this fact let us mention that the 
fine and hyperfine structure Hamiltonians in non-relativistic Quantum Mechanics are derived 
by using semi-classical arguments 12| whereas the Breit equation does so automatically in 
a systematic way without reference to semi-classical arguments. This offers more insight 
into semi-relativistic two body Hamiltonians. For instance, the standard Breit equation 
for two spin-1/2 particles can be extended to a spin-1/2 spin-0 system 13| necessary for 
exotic atoms. Higher order terms in the expansion can be taken into account and as 
we will show later, the finite size corrections can be taken into account in a straightforward 
way. Indeed, the finite size corrections to the potential have to do with the way the photon 
interacts with an extended particle, i.e., with a mild modification of the vertex for point-like 
particles (which is to say that the modified vertex will now include the form-factors). Hence, 
the modification of the vertices in the elastic amplitudes leads after a Fourier transform to 
a potential which takes into account the finite size corrections. This is a straightforward 
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generalization of the standard Breit equation where the finite size corrections are included 
in a natural way using a one-photon exchange diagram. The finite size corrections to the 
potential can then be applied to calculate the finite size corrections to the energy levels of the 
electronic or muonic hydrogen atom. The necessity to do so is the ever increasing accuracy of 
the theoretical QED corrections to the energy levels and the accurate experimental results. 

The phenomenology of the hydrogen atom cannot be disentangled from the historical 
development of Quantum Mechanics (QM). Indeed, any progress in QM, going from non- 
relativistic to relativistic QM and finally to relativistic Quantum Field Theory, was accompa- 
nied by new corrections to the energy levels of the hydrogen atom and its cousin the muonic 
hydrogen (/i^p, where a negative muon replaces the electron in hydrogen). The hydro- 
gen atom thus became a precision tool for testing predictions of Quantum Electrodynamics 
(QED) and more generally of electroweak theory, including symmetry considerations fl^ . 
The experimental techniques have reached an extraordinary precision 15| which sometimes 
surpasses the theoretically calculated values expanded in the fine structure constant a. For 
example, one of the most precisely measured quantities in physics is the hyperfine structure 
(hfs) of the ground state of hydrogen atom 161], which is known to one part in 10^^. In spite 
of this progress, there still remains one component, namely, the structure of the proton in 
the hydrogen atom which introduces uncertainties in the comparison of theories with exper- 
iment. Here, we demonstrate that the precision values calculated using QED get blurred by 
the nuclear finite size corrections (FSC). We focus in particular on these corrections in the 
hfs of hydrogen and We use the theoretical framework of the Breit equation which in 
the case of the electromagnetic (EM) form factors taken at zero momentum transfer gives 
the standard hyperfine Hamiltonian. The frequencies of the hyperfine intervals for the (IS") 
and (25*) levels in hydrogen are first evaluated using the Breit equation with and without the 
EM form factors and compared with the precise experimental values. The proton structure 
correction to the Sternheim hyperfine interval [l7] D21 = 8Ehfs{'^S) — Ehfsi^S) is evaluated 
and found to be much smaller than that obtained using the Zemach method 18|. Finally, 
the FSC for the hfs in ^^p which was used as an input for calculating the proton radius 
(rp) in a recent accurate claim of the measurement of |l9] are also evaluated. The above 
calculations show that the accurate measurements are not just limited due to the structure 
corrections and uncertainties in the proton form factors but also depend on the approach 
for the finite size corrections. 
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In the next section, we introduce the framework of the Breit equation and derive the 
potential with finite size corrections due to the structure of the proton. After having pre- 
sented the full Breit potential with form factors, in section III we focus on the part of the 
potential which corresponds to the hyperfine structure in hydrogen atom. The expressions 
for the corrections to the energy levels are evaluated using time independent perturbation 
theory. In section IV, we present the numerical results for hyperfine splitting in electronic 
(e~p) and muonic hydrogen and compare them with available data as well as results 

in literature obtained using other methods to incorporate the proton structure effects. The 
relevance of the results of the present work is summarized in the last section. 

II. THE BREIT POTENTIAL 

To start with, we briefly introduce the framework of the Breit equation. It involves an 
expansion of the elastic scattering amplitude, say, ep — >■ ep in powers of thus having 
the advantage of systematically taking into account the relativistic corrections. Starting 
with the one photon exchange Feynman diagram for a system of two point-like spin 1/2 
particles (like e~ ii^ for example), the full Hamiltonian for the e~p system considering point- 
like protons can be written as. 



where it is convenient to split the potential into several parts (to be discussed below): 

C>(Pe,Pp,r) = i:V^.(Pe,Pp,r)- (2) 
1=1 

The free Hamiltonian is expanded up to the same order in 1/c as the potential. In our case 
we have for example. 



^ = ^f) + ^f + t/(p„pp,r) 



(1) 



Pe,p PejP 





(4) 



The Fourier transform of C/(Pe, Pp, q) is the potential C/(Pe, Pp, r): 
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where in the center of mass system we can identify p 



-iV. The standard result 



for the potential in momentum space is 



25|, 
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Here = a. The meaning of the terms is as follows: 



Vi is obviously the Coulomb potential 



V2 and V3 are the Darwin terms which are also present in the Dirac equation 



20| 
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Vq, Vj are called retarded potentials 



V5 and Vg are the spin-orbit interaction terms which give rise to the fine structure; V4 
and Vs are also a spin-orbit terms, however, traditionally added to the hyperfine part 



2l| 



• ^10) ^11 are the standard terms of the hyperfine Hamiltonian. 

It is clear that the hyperfine part given above is valid only for point-like particles. For the 
electron-proton system we miss the correct magnetic moment of the proton. Indeed, for the 
hydrogen atom, the standard hyperfine potential as in text books |12| reads in natural units, 
^ = c = 1, 



hfs 



« 9p 



with Qp = 5.58. Replacing 



3(S.-r)(S.- r)-S.-S, _^ 8vr _ 3 



(18) 



^e,p " e,p/ 

and comparing with Eq. ( fTTl) . one can see that Vu in ( fT7|) gives the hyperfine potential, 
however, without the factor gp/2. The missing factor in the Breit Hamiltonian of Eq. ( IT711 
is part of the second form-factor of the proton, namely. 



[1 + 



[1 + F|(0)] = gp/2. 



(19) 



The purpose of this small exercise is to bring to the reader's attention that for point-like 
particles, -^2(0) is very small (and hence gp ^ lin the Breit potential in ( [T7|) ). as this quantity 
arises from one and more loops (anomalous magnetic moment). The hyperfine potential in 
(!T8|) already includes part of the finite size corrections (in the form of gp = 2[1 + -^2 (0)] )• 
To remedy the situation for the e~p system in the Breit equation we enlarge it by allowing 
the electromagnetic form-factors in the vertex for the proton (and for the sake of generality 
also for the electron) as follows: 



u{p') r^(p', p) u{p) = u{p') YFiiq') + 7^F2{q')a^''q. u{p). 



2m„ 



Doing so, we obtain the following potential in momentum space 



(20) 



13| 
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Amic^ 2mpC^q^ 2mempC^a^ ArriempC^q^ 



(^e (Tp \ ^ peppf 1 ^0-e.(qX pJ ^ io-e.(qX Pp) ((Te.q) (<Tp.q) 

Arriempc'^ J ^ ""^ \4:m'^c'^ 2m1c^q^ 2mf.mpc'^q^ Am^mpC^c^ 



Amf.nipC'^ ' 



[(Te.q){(Tp.q) 



\Am(.mpC^ Am(,mpC^q^ 

If one now takes all form-factors at zero momentum transfer: -Ff (0) = 1, -^2^(0) = k 
1159.6521811(7) x 10-^ Ff (0) = 1, F'^{<d) = Kp = 1.792847351(2), the potentials are: 



(21) 



fine 



for the fine structure and 

„2 



r X Pp) , , 2me (r X p ) \ 

+ 2Ke) f^(l + Ke) .CT, 



(22) 
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+ (1 + K,)(l + Kp) 



AmempC^ 



up 

^ + 3 + ^.e-p5(r) 
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(23) 



for the hyperfine part where we have included some part of the spin-orbit interaction terms. 
With = a and neglecting Ke and the term suppressed by merrip, the above fine structure 
term due to spin orbit coupling (Eq. |22]) reduces to 



book 
fine 



a 



ir X p. 



a 



2m1c^ \ 



(24) 



found in books 



12| (note that in 12|] the notation is slightly different and = 47reoa) . The 



hyperfine potential after neglecting the Kg term, the spin orbit term and the one suppressed 
by mp, reduces to Eq.( fT8|) which is also a standard result found in books [12|. If we retain 
the spin orbit term (not suppressed by m^) however, we recover 

3(Se • f)(Se • f) - Se • Sp SvT , 



Vi 



book 
hfs 



niempC^ 2 



ir X Pe 



YSe-Sp5=^(r) 



(25) 



which in the absence of any proton structure, i.e., dropping the factor gp/2 is another text 
book result as in 21|. 



We emphasize again that we have obtained this agreement after retaining the proton 
form factors at = 0, i.e., Ff (0) = 1 and ^^(O) = Kp which is not a small number. F^O) 
would be tiny if the proton would be a point-like particle. This implies that the hyperfine 
potential in books, necessarily includes parts of the finite size corrections (in the form of 
F2'(0)). It is now straightforward to generalize the result to take into account the full finite 
size corrections to the potentials by keeping the full dependence, i.e., using Ff (g^) and 



III. HYPERFINE SPLITTING IN e'p AND p ATOMS 

Having discussed the complete Breit potential with electromagnetic form factors, we shall 
now focus on the part of the potential which gives rise to the hyperfine splitting in hydrogen 
atom. We do this with the objective of evaluating the finite size corrections (FSC) to the 
hyperfine energy levels and comparing them with precision data which is available for the 
I — and I — 1 levels. As outlined in the previous section, normally terms suppressed 
by small factors arising due to large hadron masses in the denominator are neglected. In 
the calculations to follow, we shall neglect the term proportional to 1 / in the hyperfine 
potential. The factor l/m^ appears in the spin-orbit term which in any case is either small 
(for / = 1) or not relevant (for / = 0). 



A. Hyperfine potential witli form factors 



The hyperfine potential in momentum space contains the proton spin-orbit interaction 
and the proton-electron (or proton-muon in the case of spin-spin interaction which are 
responsible for the hyperfine structure. 

(o-x ■ (Tp) {(Tx ■ q) ■ (cTp ■ q) 



{2na) 



icrp(q X px 



,2„2 



mx mpC^q^ 



mx TUpC^ mx nipC^q 

[F^ (Ff (q^) + F^q^))] + 



[(i^f +i^2^)(W) +F|(q2))](26) 



(27rQ;) 



with, X — e, n for the electron or muon, F^^ — 1 and F^ — kx (the anomalous magnetic 
moment). Ff (q^) and F|'(q^) are the two electromagnetic form factors of the proton. Early 
experiments on electron-proton elastic scattering showed that the cross sections can be 
written in terms of two form factors, G%{q^) and G\^{q^), where q^ — q^ — c^- These Sachs 



form factors 2^ can be interpreted in the Breit frame to be the Fourier transforms of the 
charge (pc) and magnetization (pa/) distributions in the proton: 



'"^ ^ dr 



GW) = I Pc[r)e 



(27) 



They can be approximated fairly well by a dipole form 23|] as follows: 



(1 + qym^y 



GW) ^ ^ 



(28) 



where = 0.71 GeV^. They are related to Fi[q^) and -^^(g^) as, 

GW) = Fl{q-) + ^^Fi{q-), 
GW) = Fl{q') + m% 



(29) 



Picking up the first terms in the round brackets following F^Ff and -Ff-Ff in the Breit 
potential in fl2Tl) and putting them together, one gets (using ^ — which is the standard 
non-relativistic approximation to derive potentials from amplitudes |2^), 



i^f (q') - q 



Aira 
q2 



GM), 



which is the standard definition in q-space for the Coulomb potential due to the finite size 
of the proton. Thus, 



(1 + q2/m2)2 



1 + Kr 



q 



4m2 



q" 



(30) 



1 ArripKp 
1 + q2/m2)2 Ami + q^ 



where rup is the mass of the proton, (1 + Hp) = fip = 2.793 its magnetic moment and = 
0.71 GeV^. 

The nuclear finite size corrections (FSC) mentioned in the beginning are thus rooted in 
the two form factors Ff and . The nuclear uncertainty in general can be traced back to (i) 
different methods used to calculate the FSC (the Breit equation formalism being one such 
method which, however, is derived directly from quantum field theory) and (ii) the limited 
precision in the measurement of the two proton form factors. In the present work, we shall 
restrict ourselves to the dipole form of the form factors given above. The reason for using 
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this approximation is that (i) the potentials as well as corrections to the energy levels can 
be evaluated analytically and (ii) the objective of this work is to present the new method 
of the Breit equation for calculating the FSC and not explore the uncertainties introduced 



due to the use of different parametrizations 



24| of the form factors. 



Replacing (130 p in (12^ and taking the Fourier transform, the potential for the hyperfine 
structure in r-space is given as, 
ajjLp 



Vhfs{r) - 
where, [ix 



Ar^mxmpC^ 
= I + Kx, 

Mr) 
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(31) 
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'1 + mr) 



1 1 
m r 



6 



e-""' (3 + mr) 



and. 



2 2 

m r 



/2(r) = /i(r)-(mV/3)e- 

/3(r) = 1 - e-'""(l + mr) - ^ 

For / = and point like protons, i.e., replacing Ff = 1 and = Kp, the Fourier transform 
of the potential fl26l) leads to the standard text book potential [25| discussed in the previous 
section. Thus 



a 



AmxTnpC^ 



3(crx ■ r) 



(Tx-o-p Sir 

+ ^ CTx ■ crpO{r) 



(32) 



r' r-^ 3 

is similar to (!23ll except for the first two terms in ( l23l) corresponding to the spin-orbit 
coupling which is absent in the / = case. Since in the ground state (/ = 0), the wave 
function is spherically symmetric, the expectation value of the first two terms in (l32l) vanishes 
and the third term with the delta function contributes to the energy -E^j*"* of the / = 
hyperfine levels. However, once we include the form factors of the proton which are smooth 
functions of g^, there is no singularity or a delta function term and both terms in the curly 
bracket of (l3Ti) contribute to Ehfs for / = 0. Eq. ( 13T1) contains the standard hyperfine 
structure plus terms involving the FSC and anomalous magnetic moments of the electron 



or muon. 



B. Finite size corrections to the energy levels 

For any n and / = for example, in case of the hydrogen atom, we need to evaluate 
the expectation value, ([3cre ■ r)(crp ■ r)/i(r) — (cTe • (Tp) f2{r)]/r^) . Using the fact that 
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/ (a ■ r) (b ■ r) sin6 dO d(f) = (47r/3) (a ■ b), the energy of the hyperfine levels as evaluated 
from fl3Tl) for any n and Z = is given as, 

^3 / 2 \3 r(n + 1) 



rpn,l=0 
^hfs 



01 TfV 
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(33) 



where, a = l/{amr) is the Bohr radius with the reduced mass rrir and 2-F1 is the confluent 
hypergeometric function of the second kind. Replacing the series 



2Fi{a,b;c;z) = 1 + {ab/c)iz/l\) + [a(a + 1)6(6 + l)]/[c(c + 1)](2V2!) + 



in fl33|) and evaluating for n = 1 and n = 2, the energies of the (15) and (25*) hyperfine 
levels in hydrogen are given as, 
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(34) 
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These are the energies of the hyperfine levels as obtained from the Breit equation for the 
ep — )■ ep amplitude including the proton form factors Ff (q^), -FflQ^) ( Eq. ( l30l) ) and 
Ff = 1, F2 = Kg at the two vertices. In the absence of the proton FSC, i.e., for m 00 
in flM|) and not taking into account the correction due to the anomalous electron magnetic 



moment 
books 



i.e. — ^ 0), the energies in Eq. (!34l) go over to the expressions in standard text 
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2l| . In the next section we shall present a comparison of the corrections evaluated 



using the above equations with accurate data available for the 15* and 2S hyperfine levels. 



In 



19| . the size of the proton was deduced from a measurement of the muonic Lamb shift. 



The total predicted 2S(^^ - 2P(^^ energy difference is a sum of contributions from the fine 
and hyperfine splittings among others. The authors claimed a high accuracy on the radius 
deduced from this energy difference and concluded that the Rydberg constant has to be 
shifted by -110 kHz/c. In view of these results it becomes important to estimate the proton 
FSC to the hyperfine splitting relevant to this transition in the most sophisticated way pos- 
sible, within a given approach. Furthermore, one should take into account the uncertainties 
arising due to the use of different approaches. Hence, we also evaluate corrections for the 
2S and 2P hyperfine levels in muonic hydrogen using the Breit potential with form factors. 
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The energy El^^ evaluated using the potential in Eq. (13T]) for the triplet spin states is, 



^hfs 
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1 / 3\ 
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where, 



and 



x{f{f + + 1) --)[!- g,{ma)] 



gi{ma) = f (ma) + 2maf (ma) + 3m a f (ma) 



g2{ma) = 2m a f (ma) with f{ma) = 1/(1 +ma). 

The quantum numbers /, j and / refer to the total angular momentum of the system, 
f = j + Sp, the electron total angular momentum, j = I + Se and the orbital angular 
momentum I. The expressions for any n and I contain the hypergeometric functions 2F1 as 
before. Since these expressions are extremely lengthy, we present them in the appendix and 
here we write the result in the form of Eq. fl35l) using a truncated expansion of 2-^1 as done 
before Eqs 

IV. RESULTS AND DISCUSSIONS 



The frequencies of the I = hyperfine intervals in hydrogen atom have been measured 
quite accurately. The precise value of the (25") splitting in literature is 177556834.3(6.7) 
kHz H and that of the (IS) is 1420405.751768(1) kHz 2^. In Table I we present the 
corrections to the frequencies calculated using Eqs for the hyperfine energies. Denoting 
the hyperfine energy with Kg — )■ and m — )■ 00 as one with only m — )■ 00 as E"^"^^^ , 
the FSC correction introduced due to the proton form factors is denoted as AE^^^ = 
EJl^g — E^°^^^ . The correction due to the proton form factors as well as the anomalous 
magnetic moment of the electron is denoted as AE'^°"' = E'j^^^ — E^ . 

One can see that the FSC are large and relevant considering the accuracy of the exper- 
imental numbers. The FSC are often included via a multiplicative factor containing form 
factors 



1^ 



28 



291 or more simply by introducing a correction depending on the proton 



radius j30| . One expects [3li the difference D21 = 8 E^j^ — Elj^ to be free of proton finite 
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TABLE I: Frequencies in kHz of hyperfine splittings in hydrogen atom. 



(nl) 


f 


^fFSC 
due to proton structure 


/\ J corr 

due to Ke and proton structure 


(IS) 


1418840.09 


-37.696 


1607.665 


(2S) 


177355.01 


-4.712 


200.958 



size corrections. From Eqs (l3ll) one can see that this difference is finite but small and we 
find ADoi = -0.0833 Hz. 



A. Comparison with other methods 

The need to include structure corrections to the hyperfine energy levels due to the finite 
size of the proton in the hydrogen atom was noticed about 50 years ago. After the pioneering 



work of Zemach [18|| many other calculations followed. In what follows, we first summarize 
the main aspects of the Breit potential method and then present a comparison with other 
methods in literature. 



1. Breit potential method 

This method follows the standard method of quantum field theory to calculate a potential 
as the Fourier transform of an elastic amplitude. As mentioned in the introduction, it is 
a well documented method and widely used in different branches of physics. The finite 
size corrections (FSC) to the potential stem from replacing the point-like vertex in the 
one photon exchange amplitude (the latter without form factors gives rise to the hyperfine 
Hamiltonian in the standard Breit equation) by the standard vertex which includes form 
factors. The result thus generalizes the standard Breit equation to encode all necessary FSC 
to the potential. It should be emphasized again that the standard Breit equation gives all 
operators (fine and hyperfine) which appear in the Hamiltonian for the hydrogenic atom. It 
is therefore true to say that the fine and hyperfine operators come directly from quantum 
field theory. A perturbation of these operators by the electromagnetic form factors seems 
then to be the most natural way to include FSC. In the second step we use the unperturbed 
Coulomb wave functions in the time-independent perturbation theory to arrive at the FSC 
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to the energy levels. Thus the correction to the hyperfine energy level, AEhfs, is given by, 
AEhfs = Ehfs- Kfs^'', where 

Eh}s = J ^*{r)VHfs{r)^{r)dr. (36) 

\E'(r) is the usual hydrogen atom wave function (solution of the unperturbed Hamiltonian) 
and V/ijs(r) is the hyperfine potential with form factors (as given in Eq. (l3T]) ). E]^°j^^^ is 
evaluated as in ( l36l) . however, with the potential Vhfsij) replaced by V";f/s"*(r). Note that 
this method involves a systematic expansion in a. This is to say, if an operator being part 
of the Hamiltonian without FSC is of order a", the corresponding FSC to the potential will 
also come out to be of the same order. 



2. Zemach's original work 

The approach of Zemach as presented in [18] is very different from that of the Breit 
potential method described above. The hyperfine Hamiltonian is not derived via the Breit 
equation and hence the FSC are not obtained as perturbations to the Breit equation. Zemach 
starts directly deriving the hyperfine energy shift AEz by introducing a magnetic field H(r) 
such that, 

AEz = /ii / mr){(T2 ■ H(r))^z(r) (37) 

where part of the FSC has already been incorporated in the wave function which is a 
solution of the Schodinger equation for an electron moving in the field of the distribution 
e/e(r). fii denotes the electron magnetic moment and /e(i") is the Fourier transform of one of 
the proton's form factors -Fi(g^). The magnetic field H(r) at a given point r which is usually 
the field due to a point magnetic dipole is then modified to include the Fourier transform of 
the magnetic Sachs form factor Gm(q'^) mentioned earlier in the article. As a result, Zemach 
finally obtains, 

AEz = -'^fiifi2{(Ti-cT2) J |^^(r)|Vm(r)f/r (38) 

where /ii^2 = 61,2/2^1 2 and /m(r) is the Fourier transform of Gui^'f')- 

A comparison of Eqs f p7|) and f l38p with Eq. f l5B]) makes it obvious that the Zemach 
and Breit approaches are completely different. Whereas the Breit approach involves the 
expectation value (using unperturbed wave functions) of a hyperfine potential with form 
factors, the Zemach approach has the proton structure information embedded in the wave 
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function (in terms of the Fi{q'^) form factor) as well as the magnetic field (in the form 

of GM{q')). 

In the same work, following a perturbative formalism as given by Karplus and Klein, 
Zemach obtains another expression for AEz (Eq. (5.8) in 18|) which apart from a small 
difference in the use of the reduced mass is nothing but the momentum space representation 
of ( 138|) . Performing an explicit comparison of this expression with ( 136|) by rewriting (l36|l in 
momentum space using fl26|) it is easy to see that there arise some additional terms in the 
Breit method as compared to the Zemach approach. 

Zemach further introduced some approximations and a new distribution function /em(r) = 
/ /e(r — s)fm{s)ds to obtain the well known Zemach formula, AEz = AEq{1 — 2mia{r) em) 
where (r)^™ = frfem{r)dr. 



3. Bodwin and Yennie's correction 



In 



32| . in addition to computing the recoil corrections to the hydrogen hyperfine splitting. 



the authors provided a formalism to obtain the proton structure corrections to the hyperfine 
splitting in hydrogen. The theoretical procedure involves the evaluation of perturbation 
kernels corresponding to structure dependent photon-proton vertices. A given perturbation 
kernel does not yield a result of unique order in a but the kernels are ranked in importance 
according to the largest contribution that they produce. It was found that it was sufficient to 
include the contribution of the one loop kernel to the hfs. The six term lengthy expression for 
AEhfs was analysed and the first term after some approximations was found to reproduce 
the Zemach formula. Thus the Zemach correction within the formalism of Bodwin and 
Yennie was given as. 



AE = E, 



2amr 



1 



p2 + 7' 



2\2 



1 + K 



(39) 



where Ep is the Fermi energy and we refer the reader to [32j for details of the notation. 
The authors found that the net contribution of the remaining five terms in their expression 
was small as compared to the Zemach correction. A comparison of the results of the present 
work and corrections obtained using other methods is presented at the end of this section. 
A related approach was presented in 28| where the Zemach correction was reproduced 



in the limit of large proton mass. The order corrections to the hyperfine energy levels 
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with proton structure included were evaluated from two photon exchange amplitudes. In 
Zemach's original work both one and two photon exchanges between the electron and proton 
were considered. Zemach noted that for the two photon exchange case the interaction is of 
second order in a (i.e. one order higher than the one photon exchange interaction) and that 
the contribution of this term is very small. The Breit potential used in time-independent 
perturbation theory gives the and higher order corrections to the energy from the one 
photon exchange diagram (this potential is first order in a as in Zemach's work) . Clearly 



approaches such as in 
the FSC. 



will differ from the Breit method and lead to different results for 



^. Order corrections with one and two photon exchange 

Note that the corrections obtained in the present work are the finite size corrections to the 
hyperfine part of the Breit Hamiltonian. Once the dependent form factors are replaced by 
those at = 0, one recovers the standard hyperfine Hamiltonian for point particles. This is 
different from some of the methods mentioned above which evaluate the finite size and recoil 
corrections within the same formalism. They recover the recoil correction to the hyperfine 
structure on substituting form factors at = in the full expression (which involves the 
Zemach plus recoil corrections) 



We emphasize here that the approach in the present work is conceptually quite different 
from the ones above which are based on the Bethe Salpeter equation. There one starts by 
writing the hyperfine energy contribution (at order a^) induced by a skeleton diagram with 
two photon exchange. This consists of a divergent integral, namely, 

8Za f dk 
-Ep 



where Ep is the Fermi energy and is of order a*^. Insertion of the electromagnetic form 
factors leads to the correction (Zemach term of Bodwin and Yennie), 



, „ 8Zam ^ f dk 



.2^ 



(G.(-A:^)-l)+r%t^-l 



(40) 



1 + K 

Note that the introduction of finite size effects through the electromagnetic form factors 
Ge and Gm has not changed the order of a in the energy (i.e. to say that the order 
is due to two photon exchange and not due to the introduction of finite size effects). In 
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the Breit potential approach of the present work, the hyperfine energy calculated using the 
potential for point-like protons is of order (this is the Fermi energy). The introduction 
of the potential with form factors (which is of the same order in a as the point potential) to 
evaluate the energy using first order perturbation theory leads to order and higher order 
corrections terms from the same one photon exchange diagram. The higher orders in a arise 
due to the perturbation theory approach and not due to more photons being exchanged. 
This is evident from Eqs fl34l) where the general expression for the hyperfine energy level 
for any n and / = is written in a more lucid way by replacing a series expansion of the 
confiuent hypergeometric functions. One can see (using a = l/{amr)) that order and 
higher corrections appear in the same expression obtained from a one photon exchange 
diagram. Setting m — t- oo in Eqs flMl) in order to remove the effect of proton structure 
returns the order Fermi energy. 



5. Friar's finite size correction 

In passing, we finally note that there exist some other approaches such as the one proposed 
by Friar [29]. The finite size correction within this approach seems to be independent of 
the parametrization of the form factors. However, the hydrogen atom wave function is not 
taken to be r dependent but rather at r = 0. Besides, it also appears that though the charge 
distribution of the proton is taken into account, the magnetization distribution is omitted. 



We end this section by mentioning that the correction of -58.2 kHz in 29| evaluated to the 



IS hyperfine level is very close to -60.2 kHz obtained in [33|. It differs from the correction of 
the present work but is of the same order of magnitude. Within the Breit potential approach 
we obtain the FSC of -37.696 kHz to the IS level splitting in hydrogen. Changing the FSC of 
-60.2 kHz to our value -37.696 kHz would change the total theoretical hfs value as calculated 



in |33| (Table 19) from 1420399.3 kHz to 1420421.8 kHz. This is important considering 



the exact experimental value 1420405.751768(1) kHz and the fact that radiative corrections 



range in the order lO^^-lQ-^ kHz |33| (Tables 14-15). The FSC to D21 = 8 El% - E^J, 



is also sensitive to the method used. In [31], the correction, D2i{nucl) = —0.002 kHz as 
compared to —0.0000833 kHz of the present work. Such a small correction makes D21 almost 



free of finite size effects. The experimental value for D21 is 48.923(54) kHz 3^. Fourth order 



QED corrections are 0.018 kHz [34] and comparable to -0.002 kHz, however, bigger than our 
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TABLE II: Energies in meV of hyperfine levels and splittings in muonic hydrogen atom. 



Level 




^EFSC 
due to proton structure 


due to and proton structure 


1/2 


5.70135 


-0.0280860 


-0.02144 


2P-, In 

1/2 


1.90045 


-0.0000064 


0.00110 


1/2 


-5.70135 


0.0000345 


-0.00329 


2Pq /o 

0/ 2 


1.14027 


-0.0000077 


-0.00034 




-1.90045 


0.0000103 


0.00056 


Splitting 








2S 


22.80541 


-0.112342 


-0.085753 


2P1/2 


7.60180 


-0.000041 


0.004390 


2P3/2 


3.04072 


-0.000018 


-0.000904 



estimate of FSC From Eqs f lM|) it is obvious that the FSC to D21 begin at order and 
hence are expected to be very smalL 



B. Hyperfine splitting in muonic hydrogen 

In Table II we present the corrections to the energies calculated using Eq. f l34|) but for the 
case of muonic hydrogen {jJi^p) 2S splitting and Eq. ( 135!) for the 2Pj (j = 1/2, 3/2) splitting 
in The individual level energies are also listed. The splitting 2S is the difference 

2^{j^ - 2Sf/2° and the splittings 2P1/2 and 2P3/2 are 2Pf=2^ - 2Pf=2° and 2V{j^ - 2V{j^ 
respectively. In fact, the energy of the 2S splitting is 4 E/j/s(2S(^^). In the evaluation of 



the proton radius in [19|, the values of the hyperfine splittings were taken from |135|, where 
the FSC for the 2S were evaluated using the Zemach method and those for the 2P case were 



not taken into account. Their FSC (taken from Table II of the first reference in 35|) of 
order and sum to -0.1535 meV in contrast to the -0.11234 meV of the present work. 
This result would change the input of AEW^ = 22.8148 meV (which includes the FSC 



corrections using Zemach's method) used in 19| to AEjfpg = 22.8560 meV. Correcting the 



2P hyperfine splitting in [19|, Ae'^^fs = 3.392588 meV would change to AeJ^{j' = 3.392570 



meV. These corrections are quite relevant considering the precision taken into account in 
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19| while deducing the charge radius of the proton. 



V. SUMMARY 

To summarize, the Breit potential for hydrogen atom with the inclusion of electromagnetic 
form factors of the proton is presented. This includes the proton structure corrections to the 
standard Coulomb potential plus terms such as the Darwin term (which takes into account 
the relativistic effects), the spin-spin and spin-orbit interaction terms (corresponding to fine 
and hyperfine structure) and retarded potential terms. We focused in particular on certain 
terms involving the spin-spin and spin-orbit interaction with the aim of studying the proton 
form factor effects (or proton finite size effects) in the hyperfine splitting in electronic and 
muonic hydrogen. The finite size corrections (FSC) to the potentials originate from the 
way the photon interacts with extended objects. It is therefore possible to include such 
corrections already in the standard Breit equation where the potential is of order a. It 
is important to note that the standard hyperfine Hamiltonian for hydrogen atom follows 
naturally from the expressions for the Breit potential with form factors, once the form 
factors are replaced by their values at = 0. This Hamiltonian contains Kp(= Ff(0)), 
the anomalous magnetic moment of the proton which is large due to the fact that the 
proton is an extended object. Hence part of the finite size correction is already inherent 
in the hyperfine Hamiltonian. By including dependent form factors, we just make the 
Hamiltonian more general and complete. We presented here the FSC for three cases: IS 
and 2S splittings in hydrogen atom, Sternheim's observable D21 and 2S and 2P splittings 
in muonic hydrogen. The theoretical aspects as well as numerical values obtained in the 
present work were compared with other existing methods for FSC in literature. The present 
work aims also at showing that there exist different sophisticated and consistent approaches 
to evaluate the finite size corrections. 

Our FSC of -0.11234 meV to the 2S hyp erfine splitting in muonic hydrogen is close to 



the order ot' corrections of -0.183 meV 28| and -0.1518 meV 35| using other approaches. 
As one can see, though the three numbers are of the same order of magnitude, they are not 
equal to each other. Thus, one can say that there is a small uncertainty in the calculation 
of the finite size effects, introduced due to the differences in the approaches used for FSC. 
The present work gives the analytical expressions for the hyperfine potential and energy 
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levels including the finite size corrections due to the structure of the proton. This has been 
achieved using the dipole form of the proton form factors. Note that the inclusion of proton 
structure effects gives rise to two very different length scales in the same calculations: the 
atomic one involving the Bohr radius a (which is of the order of 10^ fm) and the nuclear 
one with oc 1/m (of the order of 1 fm). This would make a numerical evaluation of the 
integrals involved in the calculations tedious and hence an analytical approach as in the 
present work is preferable. 



APPENDIX: Hyperfine energy levels with proton structure corrections for any n 

and / / 

The hyperfine energy for any n and / 7^ is evaluated by taking the expectation value of 
the potential in Eq. flHT]) . We start by rewriting the potential as follows for convenience: 



a 



hfs 



■mf,mpC^ 



3(Se.r)(S,.f)( 1 + ^) - (Se.S,)(l + 



rriempC 



(A-1) 



where, 



hi 
hi 



11 33 

-e-'""(l + mr) e"'"'^ e"'"" 

2 6 
22 33 

-e-^Vl + mr) e"^' e"™'" 



-e-^Vl + mr) 



2 2 
m r 



;i + 2ft:p)(l + /i3) + 



1 - F) 



" e-^Ul + mr) + ^^I^e"'"'' 



1 - fc2 2 



e-'"'="(l + mA;r), 



and k = 2mp/m. Eq. flA-ll) was obtained using the identity, 



(6-,.V)(6-p.V)G'(r) 



\r or 



ldG{r)\ {&,.r){&p.r) fd^G{r) ldG{r) 



r dr 



(A-2) 



Using first order time independent perturbation theory. 



111 

hfs 



Klm{^)yhfs^nlm{r)dr, 



20 



where, 



/ 2 \3 (n-/-l)! 



1/2 



\naj 2n(n + /)! 
is the wave function of the hydrogen atom. The identity, 



e-'-/"'^ ( 5) L^LY-i(2r/na) 0) (A-3) 



n 2 



r(n + Z + 1) 



"v^' ^2(n - Z - J - 1) 

22("-'-i)(n-Z-l)! j^^ 



n 



T 2(21+1) / 4r \ 

l-j-1 Jj\r{2l+j + 2) \na)' 



allows one to write the expectation value of a function sandwiched between the hydrogen 
wave functions as. 



9 X 2Z+3 

(^) = ( iz^ 



1 



na 



J 2n22(»-'-i) 



E 

3=0 



%n-l-j-l)\ 

, n-l-j-l )j\r{2l+j + 2) 



(2j 



X 



(A-4) 



This leads to. 



Etfs U, f, (Se.Sp)) = \(-{l + Ke){l + K,) + 1) - l{l + 1) - 

+ ^) + (iw + 1) + m + 1) - ^) (d ^^){^, + '^)^ 



X 



X 



2mp r-^ 



7(/ + i)-j(j + i)-i\ 



4i(i + 1) ; 

where, the expectation values are: 

21+3 



(A-5) 



3 ^3 



f2{n r(4/ + 2j + 3) 

nay 2n22('»-'-i)r(4i + 3) 



2_Y 



yj!r(2Z + j + 2) 



r(20 (^) (-2j, 2/; 4/ + 3; 2) - r(2/) ( 
-mr{2l + 1) 



'''' )''2F,(-2j,2/;4/ + 3; ^ 



2 + mna 



2 + mna 



f-2j, 2/ + 1; 4/ + 3; — ^) - ^r(2/ + 2) 
/V 2 + mnaj 2 



2 + mna 

na \ 2'+^ / 4 
2i^i -2j,2Z + 2;4Z + 3; 



2 + mna 



_!!!lr(2Z + 3) 

2 + mna J 6 



na \ 2'+^ 
2 + mna) 



2Fi(-2j,2/ + 3;4/ + 3;— ^ ) 

V 2 + mna J 



(A-6) 



2/+3 



^^3 j 2mpr^ I \naj 2n22("-'-i)r(4Z + 3) 



E 

J=0 



'2(n - Z - J - 1)\ r(4Z + 2j + 3) 
, n-l-j-l )j\r{2l+j + 2) 



-{i + 2Kp)] r(2/) ( 



na\ 



2Z 



21 



21 



(-2,. 2M( + 3; 2, - (l . . ^ (l . 2., - ^ ] ) r(20 
(-2J. 21; 41 + 3; - (l + ^^{l + 2., - ) ,„r(2; + 1) 

TT^^Y*' 2F, (-2j, 21 + 1; 4( + 3; x-^— ) - f 1 + + (l + 2«, - 

Z + mnaJ \ 2 + mnaJ \ Zmp \ 1 — 

!f r(2i + 2) (-^) Ui, 2i + 2; 4i + 3; - i!!^^^r(2i) 

2 \2 + mnaJ \ 2 + mnaJ 2mp[l — k^Y 

na { r, ■ r,i M o 4 \ rUe KpUik f na \2'+i 



U + mfcnay' V ^' ' ^ ' 2 + m)tnaJ " 2^p{l-k^)^ + ^^2 + mkna. 
2Fi(-2i,2/ + l;4/ + 3; 



2 + mkna, 



(A-7) 



and 



e— \ = r 2 f 1 "y^^ /2(n - / - J - 1)\ r(4/ + 2j + 3) 

/ \na) 2n22(«-'-i) 1^ „ - / - j - 1 ) j\T{2l +j + 2) 

r(2/ + 3) / na ^ ( , , 4 



r(4/ + 3) V2 + mna 



- 2i^J_2j,2Z + 3;4Z + 3;— . (A-8) 

u V 2 + mnaJ 
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